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A BSTRACT 


-Skadar range and elevation angle errors which arise due to 
the refractivity varia,tions of the troposphere have heen discus- 
sed in this thesis. Eight radio sonde stations, namely, pelhi 
Allahabad, Jodhpur, Eagpur, Bombay, Calcutta, Madras and 
Trivandrum are considered. Error calculations are based on the 
known average monthly surface refractivity values of these 
places'v^^ich errors are of considerable importance for long 
range surveillance radar systems which are used for the detection 
of long range airborne targets, missiles and satellitesD 

4 numerical integration technique along with the negative 


exponential profile for refractivity variation with height has 
been used for finding the errors ."^he height and the ground range 
of the target are also fo\md out by this method ^^^he surface ref- 
ractivity, Net# is converted to the refractivity at sea level, 
for avoiding the altitude dependence of Mg and for better present 
tation, A target varying in range from 10 Kms to 600 Ems and in 
elevation angle from 0.5 to 75 is assumed for finding the ^ 
theoretical values of elevation angle and range errors 

j It is found that the elevation angle and the range errors 
of the target increase with increase in range and decrease with 
increasing elevation angle. One way maximum radar range error is 
found out to be Varying from 114 meters to 119 meters for Delhi, 
Magpur, Allahabad and Jodhpur for July /August, 116 meters for 
Madras for October and Trivandrxim for May, and from 119 met ers 
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to 120 meters for Bomtay and Calcutta for June. The elevation 
angle error is foimd out to be nearly 0.4° for Delhi and Nagptir 
for September, Mlahabad for August, Jodhp-ur for July, Madra.s 

and Bombay for April, Calcutta for May and Trivandrum for June. 

■ . . 

"''’"'Monte- Carlo technique is used for numerical simulation 
of troposphere over Delhi by generating 1000 normally distributed 
pseudorandom (IT^) numbers and these are then used for finding the 
elevation angle and the range errors for the target. Mean and 
standard deviation for elevation angle and range errors have 
been calculated and their density distributions are given, 
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DfTRCDPCTICaT 

1«1 Bac3igrotiD.d and Motivation 

An important consideration which must be taken into account 
for any radar system is the degradation in the signal-target in- 
formation due to the refractivity variations in the earth *s 
atmosphere. Signal deterioration is caused by spatial inhomo- 
geneities in the atmosphere which are continuously varying with 
time. These spatial variations produce statistical bias errors 
while the time varying components resxilts in fluctuating or rms 
inaccuracies. These effects result in refractivity bending/ time 
delays, doppler errors, rotation of the plane of polarization, 
dispersion effects and the attenuation of the signal. The first 
two of these result in elevation angle and range measurement 
errors of the target, respectively, and these two are analysed in 
this thesis. 

In general, the regions of the atmosphere which affect the 
propagation of em waves are the troposphere and the ionosphere. 
Troposphere is a nondispersive medium as the refractivity varia- 
tions and hence the tropospheric refraction effects are frequency 
independent upto 100 GHz [l], thus affecting the performance of 
most of the operational radar systems . Such effects are relatively 
unimportant for short r^ airborne radar systems such as weapon 
control radar, side looki^ terrain mapping radar and weather 
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avoidance radar# On the other hand, tropospheric refraction 
effects are of great importance in long range surveillance radar 
systems. However, on the other hand, ionosphere is a dispersive 
medium as its refractiye index is frequency dependent and at 
sufficiently low frequericies, the ionosphere is birefrigent and 
anisotropic and a simple Snell’s law approach is invalid. 

While a lot of work has been done in the other countries 
in this field, particularly by Fational Bureau of Standards, 
TJ.S.A.i hardly any work has so far been reported in India. In 
the ensuing discussion, the ionospheric effects have been ignored 
by assuming a sufficiently high frequency which is unaffected by 
ionospheric refractions. This frequency is greater than lOOOMHz 
at day time and a few hundred MHz at night time. The resulting 
elevation and range mmsurement errors are then due to the 
tropospheric effects alone and are analysed for a normal, static, 
ductless, undisturbed atmosphere, based on monthly mean of 
refract ivity for various places in India. 

1.2 Radio Refractive Index of Air 

Since radio energy at frequencies greater than 30 MHz 
is not normally reflected by ionosphere, the variation in the 
characteristics of the received fields is attributed to the 
Variations in the lower atmosphere and in particular, to the 
variation of radio refractive index of trcposphere. Electro- 
magnetic waves travel in strai^t lines in free space, however, 
in the atmosphere, there is a normal downward bending of such - 
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rays and is attributed to the normal decrease of refractive 
index, n, with height. Sadar range and elevation angle errors 
are determined by the time dejay and the bmding of the radio 
ray# respectively, which are dependent on the refractive 
index. In view of this, it is necessary to evaluate refractive 
index of the medium in which the rays are travelling. 

Refractive index can be measured directly or indirectly. 
The direct measurement is based on the relation n = “ , where- 
C is the velocity of light in vacuo, and v is. the velocity of 
em waves in air, n is the refractive index. Radio frequency 
refract ometers which are sensitive to the speed of propagfition, 
are used for direct measurement, however these are complex 
and expensive and are not commonly used. Indirect method is 
based on the meteorological variables like water vapor pressure, 
air temperature and air pressure. Such methods are much simpler 
and are mostly used all over the world. In troposphere, the 
refractivity, N, at a given height h, is given by 

H » (n-l)xlO^ = |(p + ^) (1.1) 

where, T is the air temperattxre in °K, p is the total 
atmospheric pressure in .mbars and e is the partial pressure 
of water vapour in mbars at height h and a and b are constants. 
Most commonly used and CCIR recommended values of constants a and 
b are those given by Snith and ¥eintrau3» [2], which are 
77.7°E/mbar and 4810°E respectively . These parameters are 
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available frcm weather service data and ba^ed on these, IT values 
can be found for corresponding heights. As per Smith and 
Weintraub [2], equation (1.1) is independent of frequency in the 
100 MHz to 10 GHz range and is valid upto 3® SHz with 0,5 percent 
accuracy in IT units, and normal encountered ranges of temperature, 
pressure and humidity. These parameters can also be found accura- 
tely from atmospheric emission data by using radiative transfer 
equation [27] but the procedure is costly and complex and is not 
very popular, 

1 . 3 Surface and Sea Level Refractivity 

After IT has been found for various heights, the next 
problem is to evaluate bending of the ray. Bending effects can 
be measured either by dividing the atmosphere into layers of 
varying thickness or by assigning a specific functional form to 
H, as discussed in next chapter. For these, the refractivity of 
ground Hg, at the place of measurement is essential. Eulshreshta 
and Chatter jee [ 3 ] *have computed Hg values for India based on 
the five year averages of meteorological observations for 
36 stations. Due to the dominating influence of the changes 
of density with height, the surface refractivity is dependent 
on the station altitude. This dependence has been eliminated by 
Venkiteshwaran etal [4] by reducing Ng to refractivity at 
sea level, and Fq charts have been prepared by them on monthly 
basis for whole of India. They have also enumerated the 
difference between formulae used for conversion of Hg to F^for 
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Indian peninsula and that given for U.S.A. hy Beam and 
Button [2], main difference being in the height of tropopause 
and in the definition of reference atmosphere- 

1.4 Reference Atmosphere 

Bolukhanov [5] has given International standard atmosphere 
now often called standard troposphere as defined by International 
Commission for Aeronavigation as the one having the following 
characteristics 

Sea level atmosphere pressure = 1013 mbars 
Sea level temperature =15*^0 

Relative Humidity =60 percent 

It is also assumed that standard pressure lapse rate and 
standard temperature lapse rate are 0.12 mbar/meter and 
0. 005 5 /meter of height, respectively. Relative humidity is 
assumed constant at all heights. Standard troposphere extends 
upto ii “fes* However, this is not applicable to Indian tropical 
climate. A standard atmosphere for tropics as motioned by 
Srivastava [6] has following parameters 

Mean sea sevel pressure — 1013-.2 mbars 

lapse rate in temperature in the lowest 5 Ens is taken j 

0.0054°G/meter. Tropopause for India [7] is 16 Ems during :| 

>nonsoons. Due to very large variation in relative humidity j 
from place to place in tropics, the tropical standard^^ ^ | 
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atmosphere defined above does not contain specifications of 

relative humidity and as pointed out by Yehkateshwaran etal [4]# 

the maximum error for converting to N is not more than 

bo 

units in any month over India, particularly so when station 
altitude above sea level does not exceed 500 meters, 

1.5 Elevation An^e and Radar Range Errors 

The elevation angle and range errors arise due to the 
refractive bending and the retardation of the velocity of 
propagation of em waves in atmosphere. These errors are known 
bias errors since their magnitudes are calculable. As per 
Nathanson [8], the uncertainities in their exact magnitudes, 
about 10 percent, from simplifying assumptions concerning tho 
state of the' atmosphere is partly from unknown bias errors and 
partly from rms inaccuracies. Calculations of these errors is 
generally based on averaged refiactivity data, data from radio- 
sonde stations or on functional expressions of refractivity with 
respect to height. For the first technique as well as for the 
second technique, layer approach is used which requires a lot of 
refractivity data for innumerable stations for various heights. 
This data is insufficient for India. In fmcticnal approach 
general technique employed was to assume a refractivity model 
in which refractivity decreases linearly with height, called^^ 
effective earth’s radius model. However, as has been well 
documented [2], [9], this leads to errors at long ranges and 
low elevation angles for which accurate results are greatly 
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desired# Better resiilts are obtained "by using one or two term 
exponential models based on averaged refractivity values, wMch 
obviously have limited applicability for a specific time and 
location, Majumdar flO] has discussed a biexponential model for 
inland regions and a two part model consisting of a linear vari- 
ation upto. 1,5 BnSi followed by an exponential variation for 
coastal areas for India. However, Srivastava and Pathak [ll],[l2], 
have deduced exponential models for inland and coastal stations 
of India for CGIR recommended months of February, May, August and 
November, which corresponds to Winter, Stammer, Monsoon and post- 
monsoon periods. However, these four months are unable to repre- i 
sent varying .climatical conditions of tropical climate like that 
of India and as such Venkateshwaran etal[4] have taken monthly | 
averages for various places in India for finding refractivity 
pattern. They have also found it reasonable to assume a negative 
exponential model for practical applications in India, and the same' 
is used for finding radar measurement errors. The results are 
accurate for inland stations and given fairly good values for j 

coastal stations of India. 

■ Tropospheric conditions over Delhi were simulated on the j 
computer by random generation of 1000 normally distributed Ng 
values by Montei^ Carlo techniques by using known mean and standard J 
deviation of ■^'average monthly Ng values. These are then used to J 
find the elevation angle and range errors of a target for 
comparative study. | 
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CHAPTER II 

DSTERMINAIIOH OP RAPAR MEASHREMEHT ERRORS 


2.1 Introduction 

In general the data available to a radar from a target is 
the apparent elevation angle (©^) and apparent range (r), refer 
to Figure 2.1, whereas the desired parameters are the true ele- 
vation angle (p) and the true slant range (R^) • When an em wave 
propagates through the atmosphere, the wave -undergoes .a change of 
direction, or refractive bending according to Fermat^s principle 
of least time which states that the path chosen by a ray joining 
two points is that path which can be traversed in the least 
possible time. !Hiis results in an error in the measurement of 
elevation angle of the target and is denoted as, 

c * - p (2.1) 

The range error, which is the difference between the 
apparent and time range of the target arises due to (a) the in- 
crease in the path length due to the curva-ture of the ray path | 
and (b) decreased propagation velocity in the atmosphere as 
compared to the free space velocity. In general (b) is less 
effective as compared to (a). The geometrical distance along the 
ray path is V . ^ 

: -V V r. , ^ ^ ' ' ■ v -i 

R = f dr (2.2) 

' . : r ' "Sin Q ■ ■ : 

: Q : 

Whereas, the radio range or electrical distance along -the ray 
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path measured by the radar that assumes the wave propagation 
speed to he that of light is, 
r 


ss nR = / 
^0 


n 


Sin © 


dr 


(2(5) 


Prom equations (2.2) and (2.3) we get one way radar range error 
of the target 


AR = R^ - R^ ■ (2.4> \ 

Once the errors given hy equations (2.1) and (2.4) are 
evsl.uated, true fix of a target can he obtained. However, the 
basic requirement is that of ray tracing and in the following 
sections, various techniques have been given for doing so. 

Also a number of represenrative methods for finding the errors 
are reviewed. 

2.2 Tronosnberic Refraction Effects on EM Waves 

The geometry of the ray path in normal atmosphere is given 
in Pigure 2.1, which also defines the variables of interest. 

X, the bending of the ray, is the total angular refraction of 
the ray path between two points. Refractive index, n, is slightly 
greater than \mity near the earth’s surface and approaches unity 
with increasing height. Assuming that refractive index is a 
function only of height above the surface of a smooth, spherical 
earth, i.e. it is horizontally homogeneous, the path of radio 
rays obeys Snell’s law for polar coordinates and is given by, 

n 2 ^2 cos ©2 ~ ^’'l ^1 ^ (2.5 ) 
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where subscripts ’ 1 ’ and ’ 2 * correspond to initial point and the 
final point of the ray path- x, as per Bean and Button [ 2 ] is 
given by, 

^12 = - / ® ? ( 2 -®) 

where n^ and n^ are corresponding refractive indexes at points 
1 and 2 and 9 - is the elevation angle at final point of the ray 
path. However, ray tracing is subject to the following restric- 
tions: 

(1) The refractive index should not change appreciably in a 
wave length, and 

(2) The fractional changes in the spacing between neighbouring i 
rays, which are initially parallel, must be small in a wave; 
length. 

Condition ( 1 ) is violated in case of discontinuity in refractive 
index, which is usually rare, or if the gradient of refractive 
index, dn/dr, is very large, in which case condition (2) is also | 
violated. In case of ducting both these conditions are violated, j 
Bean and Button [2] have found out that ducts do not occur for 
altitudes greater than 5 Ens and for more than 15 percent of the ] 
total time , which has also been confirmed by Chatter jee [15] for 
various places in India. In addition to these restrictions, the 
requirement of horizontal homogeneity in order to derive equation j 
(2.6) is not realized under actual atmospheric conditions. How- 
ever, Bean and Button [2] have shown that this effect can be i 
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ignored between sea level stations, which are at a distance of 
less than 500 Kins from each other and also horizontal variations 
are most effective along the surface of the earth, thus affecting 
only those rays which follow a path close to the earth^s surface. 

2.3 Evaluation of Ray Bending 

*''12 given by equation (2,6) cannot be evaluated easily 
and for solving it,, two approaches have been followed; one is the 
approximation to the integral and the second is the approximation 
to the atmosphere, fhe first leads to fimctional approach and 
the second to layer sipproach. Based on these, following methods 
are generally employed 

(1) Schulkin's method 

(2) Effective earth’s radius model 

(3) Modified effective earth’s radius model 

(4) Exponential model 

(5) Initial gradient method 

(6) Departure from normal method 

(7) Regression analysis 

(8) Graphical method. 

2 . 3.1 Schulkin’s method : 

If K-profile is known extensively, then this method gives 
accurate values with an error less than O.O4 percent. The R- 
profile consists of a series of values of R for different heights 
and a linear variation is assumed between the tabulated points. > 
The R versus height profile consists of a series of interconnected 
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linear s events, and equation (2.6) is then integrated over each 
segment, yielding the desired results. If in denominator of 
equation (2.6), n is dropped, which gives the ahovementioned 
0,04 percent error, then segmentwise integration yields 

^2 

^(rad) = - / cot © dn 

^1 

tan 9-1 + tan ©* 

1 2 


or, ^(mrad) 


2 ( 11 ^ - ^^ 2 ^^ 
tan + tan ©^ 


(2.7) 


which is a result having an accuracy much better than the desir- 
ed value, so, tan © is replaced by © for values of © less than 
10*^, resulting in an expression. 


^T^^(mrad) 


2(Ni “ ^2 ) 

61+ 02 


( 2 . 8 ) 


where now 6-^, 02 are in mrads. This approximation introduces 
an additional error nearly equal to 1 percent. The bending for 
the whole profile is then 


. T„ (mrad ) ^ I a . a ~~ — 
^ k«l ^+1 


(2.9) 


This is a very acc\jrate method and other methods try to achieve 
its accuracy. Main drawback is that a vast amount of data is 
required which is not available in India. 
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2«3»2 Effective earth’s radius model 

This method is due to Schelleng, Burrows and Eerrell and 
takes the "bending of the ray path into accoTint hy defining an 
effective earth’s radius, a = ka, where a is true earth’s radius 
and k 4/3* The. ray then assumes a straight path and various 
parameters of interest can "be calculated easily. Equation (2,6) 
after removing n, and su"bstituting dn/dh = -l/4a [2], "becomes, 

( 2 . 10 ) 

" 1 

where 6 is determined from the relation 


«h = (®o + i i (2-11) 

and for the case h^^ = h^ = G and 0 ± 10°, equation (2.10) is 

written as, 


h 



o 


dh 

4aO 


This model also gives radio horizon as, 


( 2 . 12 ) 


Goh ” VCSkah) (2.13) 

and if is in miles, a in miles and h in feets, then equation 
(2.13) reduces to 

^ ^ = y’2h (2.14) 

Here h is the antenna height. 

Extensive amoiHit of tables, graphs and maps of gradients 
of N that permit world wide estimation of the seasonal range of 
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k are available and from these radio horizon distance can be foxmd 
out, and other parameters are calculated. However, it is well 
documented that this model does not give very good representation 
of actual atmospheric H-structure above a height of 1 Km. 

2 . 3 • 5 Modified effective earth’s radius model 

The first term of equation (l.l) involving P/T comprises 
atleast 70 percent of the total value of and is proportional 
to air density, which is a well known exponential function of 
height, so one expects refractivity to decrease exponentially with 
height. Taking this into accoxont the deficiency of effective 
earth’s radius model can be reduced to some extent. Since the 
first term in the expression for refractivity (equation (l.l)) is 
proportional to air density, and the water vapor term is negli- 
gible at an altitude of 16 Em for India [7], the refractivity also 
tends to be constant at this altitude and is equal to 46H units. 
So, linear model is used upto 1 Kin and N can be represented by 
exponential function from 1 to 16 Kins, giving the expressions 


H(h) = I^fg + - iig)AN 

hgi. h< hg+1 

(2.15) 

C (h-h tD 

and H(h) = % e 

hg^^< h< 16 Kins 

(2.16) 


where AH = Hg - and its values are'hvailable for va,rious 
Indian stations, ^ 

: " ^ H. . 

Also, C = --i—- In — (2.17) 

: : : ® ; \ 15 ^^^ : ; 46 . 

where is the value of F(h) at 1 Km altitiKae, 
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Above 16 Zms, where less than 10 percent of total bending occurs, 
a single exponentia.1 decrease of U can be assumed, as given by 
the expression [7], 


= 700 


h ^ 16 Em 


{2.18,) 


The three part model expressed by equations (2.15), (2.16) and 
(2.18) has the advantage of the effective earth's radius approach 
plus being in reasonably good agreement with the average F structure 
of the atmosphere. 


2.5*4 Exponential model 

The next better approximation is the one in which a single 
distribution of H is assumed and this is called negative exponen- 
tial model, where 

-C (h - h„) 

E(h) = Eg e ® ^ (2.19) 

and =* In — § — - = In — (2.20) 

® K(lKm) Eg -A E 

Equations (2.19) and (2.20) are used to determine E at all 
heights. This model is a close representation of the average 
refractivity structure for all heights for most of the regions 
of the world . Further, this is an entirely continuous func tion 
and is extensively used in theoretical studies all over the 
world. However, data for a sufficient period of time should be 
available to find average Eg and 4 E f or better accxiracies. 
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2.3*5 Init ial gradient method 

The effect of anomalous initial N-gradients on ray propa- 
gation at near zero elevation angles, and for gradients less than 
ducting i.e., dU/dh > -157 N units/Km, may be very large. This 
may give errors if exponential model is used. To avoid this error, 
the initial gradient method utilises both the initial value, 

b 

and the initial height gradient, 61T‘, of the refractive index 
over first 100 meters above the earth's surface. Basically, it 
is assumed that exponential model is applicable for all heights, 
and most of the error is because of the difference of 6N* from 
the value predicted by the exponential model. Total bending, t, 
at any height is given by 


In practice 6N' is approximated by 6N, where. 


( 2 . 21 ) 


( 2 . 22 ) 


i.e* 


H - IT 

S % 

h - h„ 


(2.23) 


Here, value of is found either from the first level reported 
by the radio sonde above the surface or from tower or tethered 
baloon measurements, where h-hg - 30 to 300 meters. For 

further analysis, it is necessary to know what value of the 
observed value of 6lT will correspond to in the. exponential^^^^^^^ 
reference atmosphere, denoting it by Hg and using equations^^^^^^^^^^^^^^^^ 

(2.19) and (2.23 )r 
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-3 * 


IT* - N* e ® S 


h - h 


or 


S 


H 


-6H(h - h ) 




- 0.1 dlT 

lle^-^^e 


(2.24) 


(2.25) 


In practice, it is more convenj.ent to use the graphical relation 
■betti’'een Ng and 6H and find the additive correction factor which 
is the simple difference between the values of 100 meter layer 
bending as determined by Hg and jfg , i.e. 


61 ^„(6IT' 
exp 


e ) = 
o 


''’exp,0.1^^S’ 


9 ) - V ^ , (1I„, 9 ) 
o-' exp, 0.1'' S 0^ 


( 2 . 26 ) 


and using this equation, equation (2. 21) reduces to 

(2.27) 


It is claimed that this method is very useful as compared to 
exponential model for initial elevation angles of less than 
10 mrads and upto a total height of 70 Ims. This method also 
predicts any possible trapping of the ray. 

2 . 5.6 O epart-ure from normal method 

This method calculates bending by two ways, one by 
assuming exponential model .and second by the observed lT(h) pro- 
file ; and is primarily intended for pointing out the difference 
between the a ctv^ ray bending and the average ray bending that 
is predicted by the exponential H(h) profile. An A ■unit is 
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defined to study the difference between the two, and is expressed 
as, 

-C h 

A(]^h) = U(h) + UgCl - e ® ) (2.28) 

Prom this 6A(Ng, h) is evaluated, which is a measure of the depar- 
ture of ir(h) from the normal exponential profile and is expressed 
as 


6a( 1^ h) = A(Fg, h) - Ng 


-C h 
e 


= N(h) - Ng e ® (2.29) 

If equation (2,28) is differentiated and substituted in equation 
(2,6), value of , obtained is, 

0 xl 


* (ii) + I 




"S fcio ^ 

(rad.) (rad.) (rad.) (rad.) 


[AA(lTg)J X lO" 


N 


fc 


(2.50) 


where 


-C h 


A A(Ng) =A lir(h) +A [Ngd - e « ) 


-C h 

e 


= A N(h) + Ng Cg e A h 


(2.51) 


and Tu (h) is the tabulated value of for various atmospheres . 

©, and - are from the Nr, exponential atmosphere. Ihis model 
K. k+1 :: ■ O , 

has the same limitations as Schulkin's method for its applicabi- 
lity to India. 
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2.3*7 Regression analysis 

Kiis statistical linear regression technique has been 
developed by Bean, Gaboon and Thayer. Here 


a + b H, 


(2.32) 


■12 “ “s 

vjhere a and b are constants and can be found by standard -statis- 
tical linear regression methods, provided, large number of 
observed x -^2 ^or various values for many values of h and 
are available. Based on these observations, tables for a,b 
and standard errors of estimates of various heights can 

be made, and from these, bendings for other heights can be found 
out for any values. Such datas are also not available for 
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2*3*8 Graphical method 

This method is due to Weisbrod and Anderson and is very 
handy for computing refraction in troposphere. This is similar 
to Sehulkin‘s method, however, one difference being that small 
angle approximation, is not taken into consideration while arriving 
at the results. From equation (2,7), total bending t is then 
given by 

n R. - R 

T<mrad) = | — __ (2.33) 

k=o 5 00 ( tau + tan Vi) 

Terms for the d enominat or are determined from curves which are 
given by Weisbrod and Anderson [2]. However, this data is not 
available for India. 
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Methods for determining ele vation and Range errors 

There are various methods available in literature for 
finding these two errors. These errors can be calculated by ray- 
tracing methods [14], or by numerical integration [l5] . Formulae 
which are simpler and faster have also been published, but the 
necessary simplification has been obtained either by assuming an 
exponential profile with approximations [l6] or by neglecting 
path curvature [17], or both [I8]. Marini [l9] has given a pro- 
cedure to find errors in the ra,pid processing of satellite tracking 
da.ta and this method gives accurate results, provided the satellite 
facilities are available, Berkowitz [20] has given details 
about various errors and has also given integrated and stratified 
layer methods for calculating these errors. A two quatric refrac- 
tivity profile for predicting range error, particularly of 
satellites, is given by Hopfield [21], 

Besides, measuring the true range and the true elevation 
angle of the target, it is also essential in certain radar appli- 
ca.tions [ 22 1, such as air traffic control or vectoring of fighters 
for intercepting the bombers, to know the accurate hei^t of the 
target. Vogier [ 23 ] has given a method for determination of height 
of the target and error aggociated with it, Blake [l5] has given 
an exhaustive method which determines true height, ground range, 
elevation error and range error of the target. This method is 
described in next chapter and is used for finding errors for 
various places in India,^ ^^ ^^ ^ ^ 



22 


CHAPTER III 

A. .computer., program for error oalculation 

3.1 I ntroduction 

A niunerical method has "been used for error calculations 
and it is designed for automatic operation over a wide range of 
elevation angles while retaining computational precision. The 
main features of this methods are 

(1) No approximations are made in the integrand. 

(2) The step size employed in the integration is varied 
automatically by an adaptive integration subroutine that 
uses Simpson's rule. 

(3) Special computations are carried out for the grazing 
ray at low altitudes* 

Refer Eigure 2.1, starting at an initial angle the 
integration is performed to altitude h to determine both the 
ground range G and the radio range R (hereafter denoted as R), 
along the curved ray path. levins etal [24] have combined the 
integrating expressions for G and R. Thayelr, [l4] has employed 
a layer concept with a power law approximating function. How- 
ever, due to paucity of data of N(h) profile for various places 
in India, a continuous nonlinear index profile, namely, negative 
exponential profile is used here and it is felt that under such 
circumstances Blake's [l5] method is direct and simpler. 
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5 ^ A Numerical Integration Method 

This method calculates height h of the target, when 
and R are available. Once, height is found; ground range G, true 
range 1^, elevation angle and range errors can he calculated. 
Using equation (2.5) and expressing 9^ ©•» smd hy 9^, the 

expression is reduced to 


cos © = 


°0 ^0 °°^ % 
nr 


(3.1) 


From Figure 2.1, as r = r^ + h, equation (3.1) reduces to, 
n cos ©_ 

OOS © = — (5.2) 

n(l+ ~ ) 

0 

where h is the ray height above the sea level. Snell’s law leads 
to a ray tracing integral which allows computation of ray path 
length R-j_» if the ray terminal height h^j^, ani initial elevation 
angle ©^ are given, i.e, R^ = R(h 2 _,©Q). However, this function 
cannot be integrated analytically and as such numerical methods 
are employed. The inverse problem of finding h^^, given R^ and 
9 q, i.e. h^ = hfR^^, ©q), is more difficult since the ray tracing 
equation cannot be explicitly expressed in this functional form. 
Blake [15] has developed a highly precise method for finding h^ 
by using a negative exponential refractivity profile ani perfor- 
ming numerical integration of the entire ray tracing integrand. 
Refractive index for ionosphere is assumed to be \mity with the 
result that there is no bending of ray in ionosphere. ]fe,rini[l9] 
states that almost all the radio ray bending caused by the non- 
ionized atmosphere of the earth takes place within 7© Ims 
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altitude, so exponential profile is ass\imed to tie applicatle 
upto this altitude only, and after this n is taken as unity. 

This method is exact and its accuracy is limited in machine 
computation hy the number of significant figures carried by the 
computer. 


5-2. 1 Ray tracing, method 

Refer Figure 2.1, taking small increments of dR and dh 
and assuming local elevation angle as Q, 


sin © = 


dh 


or 


dR = 


dR(geometrical ) 
hdh 


ndh 

dR 


ndh 


sin “© VIl-COS^©) 

Using equation (3*2) in equation (3*3) and integrating 


(3*3) 


©„) = /' 


ndh 


(3.4) 


a n„ cos 2 

m - (- ° 


n(l + 


FT ^ ^ 


Similarly, growid range, G, measured along earth's surface at 
■sea level is 


h. 


G 


(%» %) = / 


: dh 




n 


'll;, cos Q_ 

. ■ b-- 0 


I. 2 

) r-lJ 


(3.5) 


Ta kin g spherical geometry and using law of Cosines,, 
slant range R^ and true elevation angle p are given by 


Ro * fih-l + 4ro(rQ+h^) sin(~)) . 


(3.6) 
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R 


Rin P = f +1 O 1 

Bin p t ~ + - — J 


- . f or 9 < - 

O ‘■■^0 0 ^ ^ 


and 


COS p — I — j for Q > ^ 

^ /K /! 


R, 


'0 4 


(5.7) 


(3.8) 


The negative exponential profile used is, 

-Ojh-ilg) 


H(h) = H e 

s 


(3.9) 


where h is the height w.r.t. sea level. Further expanding 
equation (3.9) 

G .h„ ~C^h 
N(h) = Ng e ® ^ e ® 


-G h 


= V ® 


(3.10) 


where Rq is the sea level refractivity of the station and is 

the decaying constant. Antenna height is ignored here as it does 

not cause any significant error in the results. N and C values 

o e 

are taken from Kulshreshta and Chatter jee [ 5 ] and RTRC monograph 
Ho. 1 .[ 4 ], 


3.2.2 Height finding method 

Next step is to find the height of the terminal point of 
ray given R and 9^, and this reqiiires inversicsi Of equation (3.4). 
Height is thus obtained by using an adaptive integration sub- 
routine that uses Simpson’s rule and the step size is varied 
automatically. Initial step computes R^^ corresponding to h^^ 
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and for "tliis, lioighi; of unrefracted ray^ h. , is given ty 

2113*21 


\a.x ^ ^ + 2 aR sin Q- )- a 


( 3 . 11 ) 


where r^ is replaced "by a, the true radius of the earth as in 
Pigure 2 . 1 . Then a value h^ ^.^ is computed, which is the height of 
the ray in an atmosphere whose gradient is linear, and of a value 
equal to that of the earth’s surface gradient of the actual 
atmosphere.. Taking effective earth’s radius a^ [2] as. 


a. 


a[ 


C a P „ 
1 - e 2 


] 


1 +P. 


(3.12) 


where p = x lO"^ and 
o ^ 


\in = ^^^e 2ag ^ ^o^ *" ^e 


(3.13) 


From equations (3.11) and (5'?13), the first trial height h^^ then 
is 


^11 “ 


h + h . 
max mm 


(3.14) 


Taking this h^^., is found from equation ( 3 , 4 ) and compared 

with desired R^.* % 1 ’ trial value of h^^ is h^^, 

where^ • 

^^ax ~ ^min^ 


a 2 


^1 - 


(3.15) 


Again, second tentative value of R^j is found from equation 

( 3 . 4 ). Iterations are^^^t^ out as per the formula^^^^^^^^^^^^^^^^ 
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h. 


In 




^^l(n-2) “ %(n-l)^ 




and range is found successively. The process is continued till 
the difference between R^^ and R^ is equal to the desired value. 
Newton’s method for comparing these two values is employed, and a 
quantity T, where 


rn ^ %n “ % 

Rn 


(3.17) 


-N 


is compared with the number 10 , where N(=6) is the number of 

significant figures to which the agreement is desired. When 
T .1 10**^, the process is terminated. The height h-^y^ is then the 
desired height at range R = R^* This value of h ^y^ is then used 
to find ground range G from equation (3.5) and R^ and p from 
equations (3«6) to (3.8). Rinally, elevation angle error is 
found from equation (2.1) and one way radar range error is found 
from equation (2.4). It is found out that for one value of 
and 20 values of range R, the computational time required is 
' 3 sec on IBM 7044. By punching 'DECK' in Column 16 of 
0IBPTC card, the whole program was converted into binary deck and 
this saved 1 mt 52 sec of compilation time of the program by the 
computer. 

As equations involved for precise numerical computation 
should avoid formulations that destroy the accuracy by taking 
the difference of two nearly equal numbers, which happens to be 
the case with equations (3*4) and (3*5) a.t small elevation angles 
and near zero heights, so these have to be manipulated into a 
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Taett^ form for accurate computation. This procedure is followed 
here. Replacing e os hy fd-sin^^S^) and assuming, 

n = 1 + X 10 vfe ® and Pq = we get from eq.uations 

(3 •4) and (3,5), the following better expressions, 


R(h^, ©o^ 


ndl+ f ) 


dh 

, f{u+V+W+W) 


(3.18) 


S(hi, 9„) = ; 


(1 + 9 q ) (cos ©^) dh 

(1 + -) X f(u+¥+¥+V.¥) 


(3.19) 


where 

-C h 

n = 1+ pQ e ® 


u = (1+ p^ ) ^ sin^ ©Q - 2 Pq 


-C.h 


V = 2p e 
o 


+ P o e 


-2C h 
o 



2h h^ 
■a a^ 


In the regions close to ©^ = 0° and h = 0, equations (3* 4) aiid 
(3.5) become infinite and as such numerical integration is im- 
possible, so some special procedure has to be adopted to over- 
come this problem. 


3,2.3 Special procedure for low angles and heights i 

For ©j^ - 0* and h = 0, the radicands in the denominator 
of the integrands in equations (3*4) and (3»5) are approximated 
by Taylor^s series upto terms of second degree, giving a general ^ 

^ ■ '4,: . ■/.''■■■'''V ■; ; . . . ; ' ■ . 4 ■ i 

: ^ :: : - , V :: : i ■ ■■ ■ : ! 
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expression of the form 

p(h) = A 

o Y.(hh^ + ch + d) 

where P represents either R or G- and A, h, C, d will have values 
corresponding to these and are denoted as d^^ and 

Aq.» bg^, and dQ_, Prom integration tables, equation (3.20) is 

written as 

(3.21) 

fb C + 2 f(bd) 

Por better computational precision, cos^^^ is taken as = l-sin*-©^ 
for Ovtieaf 0®. And the constants are expressed as, 

Ag = 1 + Po 

'’e “ ^ ■ “1”+ stn'&oClOy”- -ZCgY + 

a a 

Cj( = 2[ (2r- So - ■»■ + 5 3 

djj = dg = sin"' Sq 
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The argument of In expression in equation (5.21) is very slightly 
greater than unity for small values of h^, so, for computational 
precision, it is expressed in the form, 

^ k ■' 

p(h ) = — _[in(l + x)] 

■ iTb 

•where 


2bh^ - 2f{m) + 2 V'(b(hh| + Ch^ + d)) 

C + 2 f(bd) 


(3.25) 


The value of 71'n(l+x) is programmed for computer evaluation to 
eight significant figures by the use of the series 


ln(l + x) = X 



(3.24) 


These approximations are used for height interval from- h = 0 to 
h = 1 meter and for -0^ upto 1^. 


3.3 The Computational Method 

Ray tracing required in the height finding computation is 
done in a piecewise fashion in order to avoid Simpson’s riile 
integration of a single large height interval. The integration is 
done from h = 0 or h = 1 meter (depending on whether 9^ is or is 
not greater than 1^) to h = 10 meters in 1— meter step and from 

, V- , 

h = 10 meteisto h =100 meters in 10 meter step. The interval 
from h = 100 to h = 70,000 meters is integrated in two steps, the 
first from h = 100 to h = 5,000 meters and the second from 
h = 5,000 meters to h = 70,000 meters (or to the specified heightif 
it is less than these values). From h = 0 to 1 meter and 9^ less 
than 1^, approximations as given in para (3* 2. 3) are used. 
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Eor h. greater than 70 Kms, the assiimption n(h) = 1 is made 
and the ray follows a straight line path^ In this case, height is 
found geometrically. The radar range and ground range corresponding 
to the height increment from h^ (= 70 Kins) to the height of terminal 
point hji^ is denoted hy DR and DO. Dt is the corresponding incre- 
ment of angle 9, The local elevation angle 6 ^ (corresponding to 
©(h^)) is found from equation (2.5), 


cos = 


(1 + Pq) cos 




-C h. 
e ® ^) 


(3-25) 


and as at height = h^» n(h^) = 1, the equation (3»25) is written 
as 


cos 0^ = 


(1 + p^) cos 
— — — — 

( 1 .^) 


(3.26) 


The height h^^j is then, from the law of cosines 


h = y[ (h^ + a)^ +• (DR)^ + 2(h(, + a) (DR) sin - a 

(3.27) 

where DR = R^ - R(h.Q) 

D& = a tan(D9) = a d9, for low angles 
where, from (3.25) and law of sines. 


(1 + p.) (cos «-) (dr) 
sin D9 = - — 

^^1 + “I ) (a + ^ 


(3.28) 


All qaanti ties are known and A E. and e can he found for any range 


and elevation angle. 

The Fort rah program is given in Appendix A and has been 
used for calculation of various parameters as-given in next 
chapter for different places of India. 
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CHiPTER IT 

CALCULATIONS AM3 RESULTS 

4*1 Introduction 

Computer program developed in last chapter has been used 
to find the true range, ground range, true elevation angle, 
height, elevation angle and range errors of the target for eight 
places in India, namely. New Delhi, Allahabad, Nagpur, Jodhpur, 
Calcutta, Bombay, Madras and Trivandrum for all the months and 
the yearly mean, based on average monthly N and C data given by 
Venketashwaran etal [4]» These places are indicated in Figure 4.1 
along with their altitude above sea level and the tsrpe of climate. 
It is assumed that the input quantities which are available are 
the apparent range, R, and the apparent elevation angle, of 
the target. 

Although the computer program is suitable for any range 
from 0 to infinity and any elevation angle from 0° to 90°, the 
actual values taken for range vary from 10 Kms to 600 Ems, and 
for elevation angle from 0.5° to 75°. Minimum range is fixed due 
to computer limitations '■and max., range is fixed on the basis 
of present day radar systems available in the country. ..Limiting 
values of elevation, angle are fixed due to the limitations of the 
computer ’ -v. ; 

4.2 Tables for Radar Measurement Errors 

As indicated in Figure 4.1, the eight selected places 
are classified into three climatic groups. It was also found by 
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observation, of elevation angle and range error values that errorB 
foU.^H SimilarT general pattern and are also nearly alike for 
places under each group. Accordingly a representative location 
from each group has been selected for further analysis and these 
are, Delhi, Madras and Calcutta. Calculated values of various 
parameters for these places which are based on yearly mean values 
of Nq and C^ are given in Tables 4.1 to 4.3 for elevation angles 
of 0*.5 » 3 » 30 and 75*^. Table 4 »4 is a consolidated table 
which gives variation limits of maximum range and elevation angle 
errors along with their mean values for typicaD range of 300 Kias 

and elevation angle of 0;5°. This table also gives F variations 

■ ' " o 

and yearly mean limits. Correspondence of error values with 
values for typical climatic pattern is found to be true for 
all ranges and angles. 

4.3 Error Curves 

On analysis of the resxilts, it is also observed that 
maximum errors occur for lowest value of elevation angle, namely 
0 . 5 °, and these carry on decreasing as the elevation angle 
increases. Also, for typical value of the elevation angle, errors 
increase with increase in range and this is particularly noticeable 
for low angle values. A typical range of 3 OO Ems was selected 
corresponding to average detection range of radar systems: in the 
country and histograms of maximum range and elevation angl^were 
drawn for all the 12 months for these eight places. Figures 4.2 
to 4.9 consist of these histograms. F ^histograms based on aver- 
age monthly values of F^ are also drawn for each place. It 
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Talsle 4.1: Table of ground range, height, elevation angle and. 

range errors versus range of a target for Delhi for 
yearly mea.n Value. 


RADAR RANGE 
(MS) 

GROUND RANGE 
(KMS) 

HMIIES 

(KMS) 

RANGE ERROR 
(MRS) 

ELEVATION ERE. OR 
(DEGREES) 

ELEV = 


.50 




10 


“9.996 

0.090 

3.4 

0.010 

20 


19.992 

0.199 

6.7 

0.019 

30 


29.998 

0.317 

10.0 

0.029 

40 


39.983 

0.443 

13.3 

0.038 

50 


49.978 

0.580 

16. 6 

0.048 

60 


59.972 

0.746 

19.7 

0.057 

70 


69.967 

0.914 

22.9 

0.066 

80 


79.968 

1.094 

25.9 

0.076 

90 


89.954 

1.280 

29.0 

0.085 

100 


99.976 

1.493 

32.0 

0.094 

150 


149.968 

2.713 

45.9 

0.137 

200 


199.838 

4.263 

58.0 

0.177 

250 


249.752 

6.148 

68.4 

0.213 

300 


299.377 

8.378 

76.9 

0.247 

T50 


349.438 

10. 95 2 

83.7 

0.276 

400 


399.299 

13.906 

89.0 

0.303 

450 


449.068 

17.223 

92.9 

0.326 

500 


498.779 

20.911 

• 95.7 

0.346 

550 


548.431 

24.978 

97.7 

0.364 

600 


598.020 

29.483 

99.1 

0.380 

ELEV = 

A 

.0° 

/’ 



10 


9.957 

0.877 

3.3 

0.009 

20 


19-913 

1.767 

6.2 

0.018 

30 


29.866 

2.670 

9.0 

0.027 

40 


39.818 

3.586 

11.4 

0.035 

50 


49-766 

4.518 

13.7 

0.042 

60 


59.710 

5.477 

15.8 

0.049 

70 


69-673 

6.413 

17.7 

0.056 

80 


79-534 

7.332 

19.4 

0. 062 

. 90 


89.521 

8.366 

20. 9 

0.068 

100 


99.468 

9.368 

22.3 . 

0.073 

150 


149.027 

14.567 

27.6 

0.097 

200 


198.621 

20.137 

30.7 

0.114 

250 


248-112 

26.032 

32.4 

0.127 

500 


297.483 

32. 406 

33.4 

0.137 

350 


546.785 

39.109 

33.9 

0.145 

*1 r" Ti 

400 


395.921 

46.104 

34.1 

0*151 

450 


444 . 97 2 

53.659 

34.3 

0.156 

500 


493.935 

61. 502 

34.3 

0*159 

550 


542.722 

69-728 

34.4 

0#1d5 

600 


591.423 

78.327 

3^.4 

0*165 


Continued. 
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ELEV = 30.0° 


10. 

8.672 

5.009 

2.7 

0.007 

20 

17.292 

10 .,018 

4.3 

0.012 

30 

25.919 

15.042 

5.3 

0.016 

40 

34.537 

20.079 

5.9 

0.019 

50 

43.125 

25.127 

6.2 

0.021 

60 

51.724 

30.197 

6.5 

0.023 

no 

60.298 

35.238 

6.6 

0.024 

80 

68,859 

40.340 

6.7 

0.025 

90 

77.407 

45.437 

6.7 

0.026 

100 

85 . 961 

50.540 

6.7 

0.027 

150 

128.489 

76.240 

6.8 

0.029 

200 

170.527 

102.233 

6.8 

0.030 

250 

212.327 

128.432 

6.8 

0.031 

500 

253.777 

155.030 

6.8 

0.031 

350 

294.890 

181.845 

6.8 

0.032 

400 

3 35.671 

208.932 

6.8 

0.032 

450 

376.112 

236.235 

6.8 

0.032 

5 00 

416.220 

263.905 

6.8 

0.032 

550 

455.992 

291.721 

6.8 

0.032 

600 

495.433 

319.916 

6.8 

0.032 

ELEV = 

75.0° 




10 

2.534 

9.650 

2.2 

0.002 

20 

5.161 

19.317 

3.0 

0.003 

30 

7.750 

§8.979 

3.3 

0.004 

40 

10.292 

38.641 

3.4 

0.004 

50 

12.846 

48.385 

3.5 

0.004 

60 

15.393 

57.970 

3.5 

0.004 

70 

17.931 

67.635 

3.5 

0.004 

80 

20.463 

77.302 

3.5 

0.005 

90 

22.936 

86.970 

3.5 

0.005 

100 

25.592 

96.639 

3.5 

0.005 

150 

37.971 

144.998 

3.6 

0.005 

2 00 

50.255 

193. 381 

3.6 

0.005 

250 

62.368 

241.789 

3.6 

0.005 

300 

74.289 

290.220 

3.6. 

0.005 

350 

86.046 

338.674 

5.6 

0.005 

400 

97.624 

387.151 

3.6 

0.005 

450 

109.057 

435.650 

3.6 

0.005 

500 

120.310 

484.170 

3.6 

0 . 005 

550 

131.422 

532.711 

3.6 

0.005 

600 

142.370 

581.273 

3.6 

0. 005 
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Table 4 

.2s Table of ground range, heiglit, elevation angle 
and range errors versus range of a target for 
Madras for yearly mean N^ value. 

RADAR 

RANGE 

(KMS) 

ground 

RANGE 

(KMS) 

HMIIES 

(KMS) 

RANGE. 

ERROR 

(MRS) 

ELEVATION 

ERROR 

(DEGREES) 

_ ELEV=0.5° 





10 

9.996 

0.093 

3.7 

0. 

010 

20 

19.991 

0.199 

7.3 

0. 

020 

50 

29.987 

0,317 

10.9 

0. 

030 

40 

59.982 

0.447 

14.5 

0. 

039 

50 

49.976 

0.589 

18.0 

0. 

04 9 

60 

59.971 

0.744 

21.5 

0. 

059 

70 

69.965 

0.912 

24.9 

0. 

068 

80 

79.958 

1.091 

28.3 

0. 

078 

90 

89.951 

1.284 

31 . 6 

0. 

087 

100 

"99.945 

1.489 

34.9 

0. 

096 

150 

149.897 

2.704 

50.2 

0- 

lil 

200 

199.855 

4.244 

, 63.8 

0. 

182 

250 

249.746 

6.117 

75.4 

0. 

220 

500 

299.651 

8.532 

85.1 

0. 

255 

550 

349.485 

10.899 

92.9 

0. 

287 

400 

399.294 

13.825 

99.1 

0. 

315 

450 

499.061 

17.117 

103.8 

0. 

340 

500 

498.776 

20.780 

107.2 

0. 

361 

550 

548.435 

24.821 

109.7 

0. 

381 

600 

598.025 

29.241 

111.4 

0. 

397 

ELEV = 






10 

9.957 

0.877 

3.6 

0. 

010 

20 

19.912 

1.767 

6.8 

0. 

019 

50 

29.865 

2 . 669 

‘ 9.8 

0. 

027 

40 

59 *815 

3.585 

12.6 

0- 

036 

50 

49.765 

4.513 

15.1 

0. 

043 

60 

59.709 

5.455 

17.4 

0. 

051 

70 

69.652 

6.410 

19.5 

0. 

058 

80 

79.592 

7-379 

21.5 

0, 

064 

90 

89.529 

8.361 

23-3 

0. 

070 

100 

99.464 

9.358 

24.9 

0 ♦ 

076 

150 

149.094 

14.555 

31.0 

0# 

101 

200 

198.646 

20.116 

34.7 

0 • 

120 

250 

248.111 

26.050 

36.9 

0^ 

134 

500 

297.484 

32.362 

38.1 

0- 

145 

350 

346,758 

39.053 

38.8 

0 • 

'A'' 

154 

*t /r*i 

400 

395.926 

46.124 

39.1 

0 • 

Ibl 

450 

444.984 

53.576 

39.3 

Oa 

166 

500 

493.926 

61.406 

39.4 

0* 

111 

550 

600 

542,745 

591.438 

69.616 

78.202 

39-5 

39-5 0. 

Continued- 
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ELEV = 

30 . 0 ° 

10 

8.652 

20 

17.291 

30 

25 . 919 

40 

34.534 

50 

43.135 

60 

51.724 

70 

60.298 

80 

68.860 

90 

77.408 

100 

85.942 

150 

128.411 

200 

170.541 

250 

212.332 

300 

253.785 

350 

294.901 

400 

335.680 

450 

376.124 

500 

416.232 

550 

456.006 

600 

495.448 

ELEV = 

75 . 0 ° 

10 

2.584 

20 

5.161 

30 

7.730 

40 

10.292 

50 

12.846 

60 

15.393 

70 

17.932 

80 

20.463 

90 

22.987 

100 

25.503 

150 

37-971 

200 

50.256 

250 

62.362 

300 

74.291 

350 

86.048 

400 

97.637 

450 

109.060 

500 

120.322 

550 

131.425 

600 

142.374 


5.003 

3.0 

10.017 

4.8 

15.042 

5.9 

20.078 

6.7 

25.126 

7.1 

30.185 

7.4 

35.256 

7.5 

40.338 

7.7 

45.431 

7.7 

50.537 ' 

7.8 

76.234 

7.8 

102.216 

7.8 

128.478 

7.8 

155.018 

7.8 

181.831 

7.8 

208.915 

7.8 

236.267 

7.8 

263.882 

7 .B 

291.758 

7.8 

319.892 

7.8 


9.657 

2.4 

19.317 

3.4 

28.978 

3.8 

38.641 

4.0 

48.304 

4.0 

57.969 

4.1 

67.635 

4.1 

77.302 

4.1 

86.969 

4.1 

96.638 

4.1 

1 * 4.997 

4.1 

193.380 

4-1 

241.788 

4.1 

290.219 

4.1 

338.673 

4.1 

387.150 

4.1 

435-648 

4.1 

484 .168 

4.1 

532.710 

■ 4.1 

581.272 

4.1 


0.007 
0. 013 
0.017 
0.020 
0.023 
0.024 
0.026 
0.027 
0.028 
0.029 
0.031 
0.033 
0 . 033 
0.034 
0.034 
0.035 
0.035 
0.035 
0.035 
0.035 


0,002 

0.003 

0.004 

0.004 

0.004 

0.005 

0.005 

0.005 

0.005 

0.005 

0.005 

0.005 

0.005 

0,005 

0.005 

0.006 

0.006 

0.006 

0.006 

0.006 



Table 4. 

3 : Table of ground range, height, elevation angle 
and range errors versus range of a target for 
Calcutta for yearly mean Nq vaRue. 

RADAE 

GROEND 

HFIIIES 

RAI^GE 

ELEVATION 

RANGE 

(MS) 

RANGE 

(IMS) 

(KMS) 

ERROR 

(MRS) 

ERROR 

(DEGREES) 

ElBV = 0 

10 

9.996 

0.093 

3.6 

0.010 

20 

19.992 

0.199 

7.2 

0.020 

30 

29.987 

0.317 

10.7 

0. 030 

40 

39.982 

0.447 

14.2 

0.039 

50 

49.977 

0.589 

17.7 

0.049 

60 

59.971 

0.744 

21.1 

0.059 

70 

69.965 

0.912 

24.5 

0.068 

80 

79.959 

1.091 

27.8 

0.078 

90 

89.952 

1.284 

31.0 

0. 087 

100 

99.944 

1.489 

34 . 2 

0.096 

150 

149.898 

2.704 

49.2 

0.141 

200 

199.834 

4.245 

62.4 

0.182 

250 

249.748 ■ 

6.118 

73.7 

0.220 

300 

299.653 

8.335 

83.1 

0.255 

350 

349.485 

10.904 

90.6 

0.286 

400 

399.297 

13.832 

96.5 

0.314 

450 

449.063 

17.127 

101.0 

0.338 

500 

498.777 

20.795 

104.3 

0.360 

550 

548.434 

24.839 

106 .6 

0.379 

600 

598.026 

29.264 

108.3 

0.395 

ELEV = 

10 

5.0° 

9.957 

0.877 

3.5 

0.010 

20 

19.912 

1.767 

6.7 

0.019 

30 

29.865 

2.670 

9.6 

0.027 

40 

39.816 

3.585 

12.3 

0.036 

50 

49.764 

4.513 

14.8 

0.043 

60 

59.709 

5.455 

17.0 

0.051 

70 

69.652 

6.^10 

19.1 

0.058 

80 

79.592 

7. 37© 

21.0 

0.064 

90 

89.530 

8.362 

22.7 

0.070 

100 

99.465 

9.359 

24.2 

0.076 

150 

149.095 

14.556 

30.1 

0.101 

200 

198.646 

20.119 

33.6 

0.119 

250 

248.112 

26.056 

35.6 

0.133 

300 

297.484 

32.370 

36.7 

0.144 

350 

346.757 

39.063 

37.4 

0.152 

400 

395. 926 

46.138 

37.7 

0.159 

450 

444.983 

53.592 

37.9 

0.164 

500 

493.924 

61.426 

38.0 

0.168 

550 

542.743 

69.638 

38.0 

0.172 

600 

591.435 

78.227 

38.0 

0.175 




Continued. 
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EI,E7 = 50.0° 


10 

8.652 

5.003 

2.9 

0.007 

20 

17.292 

10,017 

4.7 

0,013 

30 

25.919 

15.042 

5.8 

0.017 

40 

34.534 

20.079 

6.4 

0.020 

50 

43.135 

25.126 

6.9 

0.022 

60 

51.724 

30.185 

7.1 

0.024 

70 

60.298 

35.256 

7.3 

0.026 

80 

68,860 

40.338 

■ 7.4 

0.027 

90 

77.408 

45. -^32 

7.4 

0.028 

100 

85.942 

50.537 

7.5 

0.028 

150 

128.410 

76.235 

7.5 

0.031 

200 

170,540 

102.217 

7.5 

0.032 

250 

212.331 

128.480 

7.5 

0.033 

300 

253.784 

155.020 

7.5 

0.033 

350 

294.899 

181.834 

7.5 

0.034 

400 

335.678 

208.919 

7.5 

0.034 

450 

376.121 

236.271 

7.5 

0.034 

500 

416.229 

263.887 

7.5 

0.034 

550 

456.003 

291.764 

7.5 

0.034 

600 

495.444 

319.897 

7.5 

0.035 

ElEV = 

75.0° 




10 

2.584 

9.677 

2.4 

0.002 

20 

5.161 

19.317 

3»3 

0.003 

30 

7.780 

28.978 

3.7 

0.004 

40 

10.292 . 

38.641 

3.8 

0.004 

50 

12.846 

48.305 

3.9 

0.004 

60 

15.393 

57.969 

3.9 

0.005 

70 

17 . 932 

67.635 

3.9 

0.005 

80 

20.468 

77.302 

3.9 

0.005 

90 

22.936 

86.965 

3.9 

0.005 

100 

25.503 

96.628 

3.9 

0.005 

150 

37.971 

144.997 

3.9 

0.005^ 

2 00 

50.256 

193.381 

3.9 

0.005 

250 

62.362 

241.788 

3.9 

0.005 

300 

74.230 

290.219 

4.0 

0.005 

350 

86.027 

338.673 

4.0 

0.005 

400 

97.636 

387.170 

4.0 

0.005 

450 

109.078 

435.649 

4.0 

0.005 

500 

120.321 

484.169 

4.0 

0.005 

550 

131.435 

532.710 

4.0 

0.005 

600 

142.373 

581.272 

4.0 

0.005 


Table 4.4: Consolidated table showing variations of R and s for various 

climatic pa-tterns. 
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FIG. 4-2 HISTOGRAMS OF fa) MAXIMUM ELEVATION ANGLE ERROR 

(b) MAXIMUM RANGE ERROR, AT RADAR RANGE OF 300 kms 
AND ELEVATION OF 0-5° AND (c)No FOR DELHI, 


No CNunts ) AR (Meters) 



FlG.4-3 HISTOGRAMS OF (a) MAXIMUM ELEVATION ANGLE ERROR 
(b) MAXIMUM RANGE ERROR, AT RADAR RANGE OF 300 k.' 


AND ELEVATION ANGLE OF 0-5° AND(c)No FOR ALLAHABA: 


No (N unts AR (Meters)- — — £ (Degrees) 
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FIG. 4-4 HISTOGRAMS OF (a) MAXIMUM ELEVATION ANGLE ERROR 

(b) MAXIMUM RANGE ERROR, AT RADAR RANGE OF 300 kms. 


AND ELEVATION ANGLE OF 0-5“ AN D (c) No FOR NAGPUR. 


No(Nunts) — — £/r3g 





FIG. 4-5 HISTOGRAMS OF(a)MAXIMUM ELEVATtOH AliGLE ERROR 


(b) MAXIMUM RANGE ERROR, AT RADAR RANGE OF 300 kms. 
AND ELEVATION ANGLE OF 0-5^AND(c) No FOR JODHPUR. 


y n ts) — — ► A R ( M e t e r s ) 



FIG. 4-6 HISTOGRAMS OF (n) MAXIMUM ELEVATION ERROR (b) MA,XIMIJM 


RANGE ERROR, AT RADAR RANGE OF TOOkms.AND ELEVATION 
ANGLE OF 0 . 5 ° AND (c)No FOR CALCUTTA. 
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FIG, 4.7 HISTOGRAMS OF (a) MAXIMUM ELEVATION ERROR (b) MAXIMUM 
RANGE ERROR, AT RADAR RANGE OF 300 kms.AND ELEVATION 
ANGLE OF 0.5° AND (c)No. FOR BOMBAY 
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FIG.4-8 HISTOGRMS OF (a) MAXIMUM ELEVATION ERROR (b) MAX- 
IMUM RANGE ERROR, AT RADAR RANGE OF 300kms,AND 
ELEVATION ANGLE OF 0-5“ AND (c) Nq FOR MADRAS 
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F[G.4-9 HISTOG^^^^ MAXIMUM ELEVATION ERROR 

(b) MAX MUM RANGE ERROR, AT RADAR RANGE OF 300 

trivandrJm"'™* M". ™r“ 


■ 50 

was also observed, that general shape of the. histograms remains 
the same for all other va-lues of radar range. 

A computer program was developed to plot the. errors versus 
radar range curves on IBM 1800. Maximum and minimum errors a^long 
with the errors based on yearly avera,ge of were drawn for all 
the eight places for ranges varying from 10 Em to 600; Em and for 
elevation angles of 0.5*^? 15° a.nd 75°. These curves were smoothe-ned 
wherever required, particularly at low error values of elevation- 
angle error curves so as to avoid the inherent incapability of 
IBM 1800, These curves are' given from Figures 4.10 to 4.17. 

4 . 4: Radar Gov e ra,ge D iag rams • 

Bange-height-angle charts are important for many radar 
applications, Eadar coverage dia.grams are plots of radar range 
versus elevation angle, and form a part of such general plots. 

A version of these charts is drawn in Figure 4.18(a) for Eew 
Delhi based on yearly mean value of FT^. Slight variation in this 
diagram is; expected for individual months, this being due to 
variation of over the year. Height versus ground range curves 
for various values of elev8,tion angle are drawn in Figure 4.18(b). 

Curves for Madras and Calcutta can be similarly drawn. 

4.5 Height ve rsus Error Curves 

If the target height and elevation angle are known, then 
a quick procedure for correcting the errors can be formulated by 
drawing height versus error curves for various elevation angles. 
These curves are given in Figure 4*19 for Delhi for maximum eleva- 
tion and range . errors. Normal aircraft ceiling heights are . taken 
in these curves. Similar curves can be di.rawn for Madras and 


Calcutta 
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F!.3.4-11 (a) RANGE ERROR (b) ELEVATION ANGLE ERROR CURVES 
FOR ALLAHABAD, 
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F,G, 4-14(0 RANGE ERROR (b) ELEVATION ANGLE ERROR CURVES 
FOR CALCUTTA. ' 
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Legend ;-1. Yearly mean 2 . M 

FIG. 4-15 (q) range error (b 
CURVES FOR BOMBA' 
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Legend ; - 1. Yearly mean 2/Minfmum crrmr 3, Maximum 
FIG.4-17. (a) RANGE ERROR (b). ELEVATION ANGLE ERROR CURVES 
FOR TRIVANDRUM. 
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FI6 4-13 RADARl COVERAGE DlAGRA>^S FOR YEARLY 
MEAN VALUES OF No FOR DElHiL 
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FIG. 4-19 (a) HEIGHT Vs, MAXIMUM RANGE ERROR & 

(b) HEIGHT Vs. MAXIMUM ELEVATION ANGLE ERROR 
CURVES FOR DELHI. 
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CHAPTER 7 

HTJMERICAL SIMJLATIOH OP TROPOSPHERE OVER DELHI 
BY MOHTE-CARLO TECHHIQUES 


5 »1 In tro d ucti on 

In the last chapter the monthly a.verage values have 
hcc-n used for finding various parameters for the target . Actually 
there is a marked vari''tion in the values Td-th time and 
enough of data is not available for analysing the problem 
exhaustively. Accordingly a large number of values of Eg, the 
surface refra-ctivity, have been generated for simulation of 
troposphere over Delhi, and this is based on the law of large 
numbers and the central limit theorem which states tha,t the com- 


bined distribution for a large number of random variables tends 
to abnormal distribution and the variance is small for large 
numbers. To start with, the distribution- is assumed to be normal 
and the mean (= 332 F units ) and the standard devia.tion (= 25 N 
units) of knoTO Eg values are used for random generation of 


1000 H„ values bv Monte-Carlo techniques for the simulation of 
the tropospheric conditions over Delhi. The term Monte-Carlo is 
frequently reserved for the use of stochastic techniques for 
solving an essentially deterministic problem, and due to this 
the random numbers generated are ca,lled pseudo random numbers. 
Program, of appendix 'B' [25] is used for generating these pseudo 

random Eg values. First uniformly distributed random numbers are 
genera.ted in the interval (0,1) and these are then converted to 
numbers in normal distribution. 
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Sg'fc 'U -P for S imulation of Troposph ere 

Once Ug values are generated, tlie decaying constant per 
Em can be found' out for each of these, from the relation 

No 

C. = ln( S — -) (5,1) 

Ng - AN 

where AN is the average of the monthly vain® for COIR recomjnen- 
ded months of February, May, August and November for 00 GMT and 
12 GMT [26] . 

The computer progra.m of Appendix 'A' was modified to take 
these changes into account and elevation angle and range errors 

were foimd for different values of N„ for a target having a range 

■ b . ■ 

of 300 Ems and an elevation angle of 0.5°. This program xv'as also 
used 'to find the means and standard deviations of Ng, elevation 
angle and range errors. 

5.3 Results 

Density distribution curves for the resulting error values are . 
dra^wn in Figure 5.1, and on calculation it was found that, as 
expected, _ values follow a norm.al distribution. Means and 

siandard deviations of Ng,AR and e for these 1000 generated 
values of Ng are 

Mean value of Ng = 337 N units 

Ar = 67.9 meters 
e = 0.268° 
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Standard deviation of 


= 27 N imits 

O 

Ar = 1.3 meters 
e = 0. O 63 degrees 


Mean and standard deviation of llg are comparalDle to th.e 
known yearly a.verage values for Delhi. Mean value of e is sligh- 
tly higher than a,nd that of AR is slightly lower than the 
knoi-m yearly a,vera.ge values for Delhi and a significa,nt contri“ 
hution to this difference is due to the approximations used in 
the determination of AN. The standard devi'^tions of errors 
from, their meajis are comparable to the reported values for " 
various pla.ces in the world. 
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CHAPTER VI 

- PISCnSSIOH OE RESULTS AMP COWCLUSIOH 

Radar range and elevation angle errors which arise due 
to the refractiyity variations of the troposphere are discussed 
in this thesis. These errors are of great importance for long 
range surveillance radar systems. Errors have been calculated 
for eight radio sonde stations for all the twelve months and 
their yearly average for known monthly average values, for a 
ta,rget varying in ranges from 10 Bns to 600 Eins and in elevation- 
angle from 0,5^ to 75°* The various results of the calculations 
are given In Chapter IV. It is observed that in all cases the 
elevation angle and range errors of the target are having a 
maximum va,lue for lowest value of target elevation angle of 
0,5*^ for all ranges. These errors carrjT' on decreasing with 
increasing elevation angles of the target for ans?- typical 
target range. Also, for a fixed elevation angle of the target, 
the errors increase with increase in range. Vanious histograms 
of Chapter IV indicate the average monthly variation of the 
maximum errors for typical target range of 300 Kms and elevation 
angle of 0.5° for all the eight places. General pa,ttern for 
individual places was found out to be the same for different 
valu es of range and elevation angle of the target. A version 
of re.dar coverage diagram is given in Cha-pter IVj which can be 
used to find any one of the three par a,mG ter s provided the other 
two are given out of h, R, or h, G, Also a method of : 
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error correction is given whicli can be used for quick deter- 
mination of the true target fix in case the target elevation 
an^e a.ud the height is known. 

Monte-Carlo techniques were used for numerical simulation 

of troposphere over Delhi by generating 1000 normally distributee 

pseud orardom (l^g) numbers by assuming a normal distribution for 

known monthly average Ng values over Delhi, The results are 

given in Chapter V and the Galc^iLated means and standard devia^ 

tions of HR and e are found to be comparable to published 
o 

data and the known values for Delhi, 

In this thesis, eleva.tion angle and range errors of a 
target have been analysed for eight stations in India* This can 
be extended to cover up more regions of Indian peninsula- Other 
effects, like, doppler effect and attenuation of the signal, 
which also are caused by the refractivity variaticjn of the 
troposphere can be calculated for various places for various 
values of range and elevation angle of the. target. 



APPEliDIX 'A' 


67 


A COMPUTER PROGRAM FOR THE EETEPJ4I1UTION OF HEIGHT, 
GROUITD RANGE, TRUE RANGE , TRUE ELEVATION ANGLE, 
ELEVATION ANGLE ERROR AND RANGE ERROR OP A TARGET 
OP VARYING lUTTGES AND APPARENT ELEVATION ANGLES. 


Ig^IBETO main 

0 TROP OS PI-ERIC PROPAGATION ERRORS IN RilDAR ESiHUREITENT IN INDIA 

C CALCULATION OP ELEVATION AND HCN® ERRCES VERSUS RilNGB POR - 

C V.VRIOUS ELEVTiTION iR-JGIES OB' E.U PaTH 

0 THE RAY HEIGHT IS COMPUTED AT SIEOIPIED R^'iNGS TiElD ELEV.'EION 
0 .INGLE 

0 PCE SRiNDCRD NEG.iTIVE EXPONENTIAL PEPRACTIVITY atmosphere 

0 TIE R.1NGE iRID HEIGHT AES GIVEN IN M AND ELBVTHION SINGLE IN 

0 DEGREES 

C EBP = NS * C 10 ** -6), GR.VD = CE PER METER 

G NS .WD GE VALUES 7RE T NOEN PRai ETRC MONOGRAM NO.l BY 

C VENKITESHWAJR.1N ET.VL (PEB 1970) 

§ GB DENOTES THE REPRiCTIVITY GRTJ)3ENT PER METER 
0 NS IS THE RBPRACTIVITY AT THE SURE ICE IN N UNITS 

G Nn IS THE VAIUE CP NS COIUHIRTED TO SEA lEVEl IN N UNITS 

0 EIBVaTION is the initial ELEVi^ I on ANGLE OP RAY, 0^, IN DE-G-cBjiJS ■ 

0 RNGl IS THE TERM IRIL R.. INGE OP RAY, R, IN KMS 

0 RNG DENOTES INTEGR.lTED R.1NGE,R, IN KMS 

C XI DENOTES R\Y HEIGHT, R, IN METERS 
0 mniES DENOTES RAY HEIGHT, h, IN KMS 
0 GRG DENOTES GROUND RANGE , G, IN KMS 

C RNGERR DENOTES THE DTEPERENCB OP RAY P:1TH *ND STRAIGHT LINE 
0 DISTANCE OP THE T.iRGET , AR, IN SETERS 

0 ELEVEE DENOTES THE DIPPERfflCE OF APPARENT ELEVTiTION ANGLE /sPID 
0 TRUE ilNGIE (CP THE TARGET , e , IN DEGREES 

0 N 01 A AND NOIB DENOOE THE OUTPUT PAR THE TEES OP THE 3 RIPS ON » S 

0 RULE NUl/ERICAL INTEGRATION SUBROUTINE, S MOON . 

0 R.\DAR ilNTENNA HEIGHT IS IGNORED THESE GxiLCULAiTIONS 

0 TalR GET HEIGHT IS TAKEN WRT SEA LEVEL ■ 

G GROUND R.NGE IS DIST.NGE ALONG THE ASSUMED SPHERICAL E-GRTH TH 

C SEA LEVEL 

GCIMON /B/ HEP, RilD, GR;L), CONST , H, A 
EXTERNAL PI, P2 

H3(H1, R2, SP ) = SQRTCRT*EI+R2*R2+2i*R1*E2*SP^ 

THT { P1,-EH, GR, hi, ED ) = iROOS( PI / ( ( ■ 1 ,+RH * (EXP ( -GE * 
1HT + Hi /fo) ))) 

ilEG(A1 , B1 , 01 , IS) = (2.*Al*X-.2 ,*S(iRT(Al*C1 ) + 2 .* S(^(a1* (a 1* 
1X*X+ B1 * Z + C1))}/(B1+ 2.0* SOKTC Al* 01)) . . 

data HilD / 6 .373E+06 / 

PRINT 600 
■ xo = 100.0 

XA = 5000.0 
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ZB = 70000.0 
COBST = 0.001 
PIE = 3.14159265 
APPA = 180.0/PIE 
P/J>A = 1 .O/APPA 
H = 1 .0 
BSIG =: 6 
IL'iZF =. 10 

PEEC =■ 10.0 ** ( - BSIG ) 

1 EjlliB 500, EEP, GR/sB, MOMI , M01WH2 
n? ( REP. LI . 0.0 ) 00 10 499 
500 PORILIT (P9.6, PI 0.7, 5Z, 2A5 ) 

. PRIFI 17, MOFTHII, M0HIH2 , REP, GR:J) 

17 POIMAI ( 1 HO, 30Z, 2A5//30X, * REP = * , P9 .6 // 30X,*GR.-I)=* , 
1PI0.7 ) ■ • 

AB = 1 .O+REP 
RI = 1 .0/R.iJ) 

RS = RI*RI 
QrM = REP*GIAUD/7iB 
GS = GiiM*G.a,I 
I = R.H)^GAI4 
PA0=:1 .0/(1 .0-1) 

PAJ)1=ILID^PAC 
ELE7 = 0.5 
5 RFG1 = 10.0 
PRIITI 993, EEEV 
REN = EIEV * P.VPA 
OT = 0 

BE 101 K = 1,20 
m = OT +' 1 
Cl = 0 03 ( REF ) 

SF = 3 BT ( REF ) 

3 = SF SF 

U = .IB *AB * 3 - 2 .0 REP - REP REP 
A = AB * COE 
RFG3 = EFG1 / COFST 
mmT = H3 ( RAEl, EFG3 , SF ) 

mLlZ . = (H3 ( RAE, EFG3 , SF ) ) * AIB 

O' ( EIEV - 1 .0 ) 56 , 56 , 57 

56 Q= 2.0 * ((1.0 + S ) * RI - GBI ) 

QQ = 2 .0 * ( ( 2 .0 * Gill - RI ) * S - GMt + RI ) 

P ='GR;E * Gi'AI + GS - 8.0 * G.'iM * RI + ( 5. + S ) * RS, 

PP = GRiiE * Gi'JI - 7 . * GS + 8 . * GA'I *RI - 3 . * RS + S' * (1 0 . 

1 * GS - 8 . QM * RI - 2 .0 *GEi1E * GR/I + 3 .0 * RS ) 

IRG1 = .1RG( PP, QQ, 3, H ) 

/iRG2 = liRG (P, a, SV H) • 

GilLL L0G1PL ( ARGl, 12 , Sm ) 

OulLL LOGIPE (.RG2, 12, SBlvII ) 

EFG = COFSI * .IB * SMI / SQRI (PP) 

GRFG = OOFSP.* CP.* SD1I2 / SCRP .(P) 

Mr: 2 

GO PO 19 
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57 M = 1 ■ 

mm = 0.0 

GOTO =0.0 
19 no 29 J = 1’,2 

in ( J .EQ. 2 ) M = 2 

no 29 I = M , 1 0 

XI = ( I - 1 ) ^ 10 ( J - 1 ) 

X2 = I 10 ( J - 1 ) 

C SUBEOUIIITS SIMPS on SUIE IMJJtEEIC/ili BlTEGiL.T lOB 

c;iii Silicon (xi , X2, prbc, 15 , eihc, noi , a, ni) 

CALL S IMG on (XI , X2 , PilGC , 15, GHIHC, n02 , R, E2) 
ilFG = illTG + HinC 

29 ghhct = oano + Game 
aiiGio = aiiG 
GailGIO = GLFG 

XI = (HlLOC + HMIH ) / 2 .0 
nELTA = (m.(uix - Hf,OT ) / 2.0 
1 = 1 
LII'I = 15 
■TEST = PEBO 

IE (HMIl - XA ) 21 , 21 , 22 

21 PAIGE = 0.0 
lOIB = 0.0 

210 C^OiL SIMGOI ( XO, XI, TEST, LBI, BIGA, lOIA, HA, El) 

GO TO 23 

22 IE (HMIl - XB ) 24, 24, 25 

24 C .ILL 3 BIG 01 ( XO, XA, TEST, LIM, EIGA, lOIA, HA, El)’ 

240 C.ILL SBiCOl ( XA, XI, TEST, LB£, EIGB, lOIB, HB, El. ) 

GO TO 23 

25 gall SMC 01 ( XO, JA, TEST, LB.I, aiGA, lOlA, HA, El ) 

C.VDL SMOOT ( XA, Xb, TEST, LBI, aiGB, lOIB, PtB.^ El) 

BBLaiG = ( aiGl - RIGB - PA1GA - aiGl O ) / OOIST 

xa = XB + HAD' ■ ■ 

Eni3 = tht( a, asE, c-EaVD, XB, a,.in ) 

S13 = sil ( aD13 ) ■ 

XI = h 3 ( xa, nELaiG, SH3 ) + xb 

HI = XB / RAD 

PHic = agsii( a * niiaiG / ( ( i . o + hi ) * ( rij) + xi ) ) ) 
1 = 1 

aiG = aiGi 
GO TO 391 

23 aiG = aiGio + aiGA + aiGB 
IE ( 1 .EQ. IL1X1 ) GO TO 30 
TESTo = ( aiG - aiGi ) / aiGi 
T1ST1 ' = iiBS (TESTO ) 

IE (TlSTl - TEST ) 30, 30, 50 

30 IP ( mm - Zii ) 38, 38, 39 

38 call S IMG 01 (XO, XI, TEST, LBI, GHIGA, lOIGA, RGA, E2) 
GEIGB =0.0 

GllGO = 0.0 
GO TO 37 

39 CiiLL SIMGQl ( XO, XA, TEST, LIM, GBIGA, lOIGA, EGA, E2) 
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0.1LL SDIGOl ( XA, Zl, TEST, LBI, OHM, ITOIGB, HOB, I'2 ) 

GPlFGC = 0.0 
GO Q?0 37 

391 SIMGOIT (XOjXAjTEST, LUI, GRBGA, BOIGA, EGA, B2 ) 

G.-JLl SEIGOF (XA, XB, TEST, LD:f, GRBGB, ITOIGB, RGB, BS ) 

GEl'TGG = PHIC * E/J) * COITST 
37 GtEBG = GRIJG10 + GEl'fGA +GEITGB + GEITGO 

THEDTA = AEPA * THT(A, EBB, GS.II), XI, E/O) ) 

.ELT = XI - ‘ ■ 

TEG = SBT ( GEEG / ( 2,0 * EJJ) * GONST )) 

BBS = SORT ( alt * AIT + 4 .0 * HAD * (E.AD + ALT ) * TEG * TEG ) 
EBGEEE = REG / G GIST - EBS ' 

IB ( ELEV - 45.0 ) 360, 360, 370 
360 T-EUELV = - APPA .iESIE ( HBS / ( 2,0 * IL’J) ) - 

IZI / PiBS - ( XI * XI ) / (2.0 E;J) * RBS )) ■ 

GO TO 380 

370 TEUBLV = 90.0 - APPA * ( ;ulSIE.((( E^il) + XI ) * ( SIH ( 

IGEHG / ( R;J) * GOEST )))) / EBS )) 

380 ELBTBR = EIEV - TEUEL? 

IBIILBS = XI GOEST 
GEdG = GRITG 

PEIITT 31 , M, REGl, GEG, fflIIIES, EIIGEEE, ELBVER 
600 BailLiT (IH 25X, * TEOPOSPHEEIG PE®;iG .BIOE EiSiOPB IE E.J),Ei MS/iS, 
lUEElBETS OB’ ELS V. IT I Oil AJCD E.EIGE OB T.UiGET * / 50X, * lOEOSPHEi-S 
2ASSin,IEB 1T01EEBRA0TI7E * // . 

3 60X, LATA BOE ) 

993 .BaRl,IAT- (IHO, '7X, * EIEV = *, B5 .1// 8X, * SL BO *, §X, * R^IDAiE EA 
lEGE * ,.12X, * GEOUIIDIWGE 12X, * HMILB3 lOZ, * E.AJGE ERROR 
2*. 5X, « BIEVADIOE ERROR ^ / 22X, * (KIIS)* , 22X,* (KI.I3 -)* ,1 8X,* (Zm 
3S)* , 1§X, *(iaES)* ,14Z, *(LE&SESS)^ ) 

31 BORIIAT (1H ,10Xj 12, lOX, B6.1, 15Z, B11.3, 15X, B8.3, 13Z,B6.1 
1 , 15Z, B7.3 ) 

GO TO 555 

50 IB ( IT .GT, 1 ) GO TO 60 
XU = XI 
EEGE A= RIG 
H ^ . +1 

LEtLTA = 0.5 * LELTA 
IB ( TESTO ) 3, 3, 15 
3 XI = ZI + LELTA 
GO TO 6 R 
1 5 XI = XI - LELTA 

6 IP ( HMIN - XA ) 210, 210, 240 • . 

60 XM = XI • : ■ - - \ - 

XI = (ElTGl-ENGE) * ( XI XR ) / (EFG ~ REGIT ) + HT 
XF, = XBi. : 

HFGF = RITG 
N = n+ 1 

GO TO 6 : 

555 IB (■ ElGl ,GT. 95.0) GO TO 333 
EFGl = RIGl + 10,0 
GO T O 101 
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333 re* ( M&l .GT. 595.0) GO TO 101 
riTGl = ElTGl + 50.0 
101 coiH)iiroE 

IE (EIEV .GT. 4.9 ) GO TO 225 
EIEV = EIEV + 0.5 
GO TO 5 

225 12' (EIEV .GT. 29.0 ) GO TO 230 
EIEV = EIEV +5.0 
• GO TO 5 

230 IE ( EIEV .GT. 74.0 ) GOTO! 
EIEV = EIEV + 15.0 
■ GO TO 5 

499 STOP 
EEE 


^IBETG SUBl 

EmiCTIOB El (X) 

C BTEIGBAIil) -EOIl SIHC QT SUEHOUTBiE 

COIEOT /B/ EEE, 'RED, GE.ID, COjvST, U, A 
BB = EEE * EXP(- GPAJ) X) 

CC = X / Eu-il) 

V = 2 .0 BB + BB * BB 


W =2.0 * CC + CC * CC 
EX = SQBT (U + V + ?/ + V 

21-.= .0PBST * ( 1.0 + V ) 

&SJURK 


^lEPTC SUB2 

EUHGTIOB E2(X) 


W ) 

( 1 .0 + CC ) / EX 


G IETEGE.iBI) EOE SBIGON SUBEOUTIEE 


GOIMOIT /!/ EEE, iL’JD, GiLlB, GOHST, U, A 
BB = BEE * EXP ( - GR.ID * X ) 

GG = X / R.J) 

V = 2 .0 * BB + BB * BB 


W = 2 ,0 * CC + GO * CG 

EX = SORT ( G + V + ¥/ + V * ¥4 ) 

E2 - CCaiST » A / (U.o + GG ) * EX ) 
EETGRR 


ERI) 


^IBETG SmCOR 

SUBROUT IKE S BIG ai( Xl/XEliB, TEST A i ni , ABEA, ¥01, R, E ) 
¥01 = O'- 
.. El =10.0 
ODD = 0.0 

jm = 1 
Y ^1,0 : 

EVER = 0.0' 

AREAl = 0.0 

19 E¥BS = E (Xl)- + E (XB¥B) 

2 H = ( XE¥D « X1 ) / 

ODD = EVER +ODI) 

X = XI + H/ 2 . 
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BVBl =0.0 
D-0 3 1 = 1, Jl® 

21 = E7BH + F(Z) 

X = X + H 

3 COHTBKIE 

31 AEEX = (Em)S + 4.0 * E^/EE + 2 .0 * ODD ) * H / 6 .0 
EOI = HOI + 1 

34 R = JiBS. ((iREAl - .IRE A ) / /JiEA ) 

IE (EOI - 111,1) 341, 35, 35 

341 IE ( R - TESI ) 35, 35, 4 

35 RBIRRE 

4 iiREAl = AREA 
46 JIE? ='2 * JET 

Y = 2,0 * Y 
GO T 0 2 
EEE 

^IBETO lOGlPI 

SIBROUTIRE LOGIPL ( AGT, N3E, SUM ) 

ARUM = ■ AGT 
DENCM = 1 .0 
SUB;T = 0.0 

1 TERM = AETJI^ / EEROM 
SUIi = SUM + TERM 

TEST = iVBS (TERJI / SUM ) 

IE (TEST _ ( 10. ** ( - USE ))) 2, 2, 3 

2 RETUTiR 

3 iUMUM = -.IRIRil * AGT 
EEROM = EEROM + 1 .0 
GO TO 1 

ERE 
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.■iP2ElCDIX - B 


mOCrRMK PCS SiWBai GBJlEE.iTIOl'T OP PSEUDO E.WBOM EUIffiEES POE 


SIMUE..EIOE OP TEOPOSPHBEE OVEE BBLHI. 


^IBPTC ISilE 

0 THIS PEOGETII USES THE SUBEOUTIHE EAJ5M TO GB HE RITE 500 ?.‘JES 
0 OP PSEU-DOEAHDOM TOIIBEES IH FOEFlAL DISTEIBTJTIOH POE A SET OP 
0 VAIIiES HLIVIHG A DIE .,\H OP 332 H UHITS AHD ST.IHDAPA) DEVIATIOH 

C OP 25 H UNITS, BASED ON KNOUH Nq HONTHLY AvCEEiGE YAIUES POE 

0 DELHI 

INTEGBE TAU , ALPA,G, ENl 
HEAD I.TU 

ALFA = 2**18 +3 

0 ENl = 2**35 - 3 CaiPUTED AS BE LOU TO AVOID INTEGEE OVEEPLO?/ 
END = 2* (2** 34-2) +1 
SIGMA = 25 .0 
m = 352 .0 
N = 1 
G = 0 
TAU = 3 

1 DO 2 I = 1,500 

GAIL EilDM (SIGMA, IIU, TAU, G,N,.mPA, ENl, SN, SNl) : 

0 THIS PE0GE.ffl IS PGR PRINT EJG THE NUI,£EEES 
PRINT 100, SN, SNl 
loo POElvIAT (2E16.8 ) 

2 CaWINUE 

IF (TAU-3) 5,6,5 

5 TAU = TAU +1 
GO TO 1 

6 STOP ' 

END 

^IBPTC REDM • • \ • ■ ■ 

SUBROUT II® EADM ( S IGIvLl , MU , T.VU , G , N , .-IIP A, ENl , SN , SNl) 

INTEGER, TAU, AIPA, G, EN, ENl 

EEAL'MU 

C NEXT - THEBE STATEMENTS PEaDUOE INTEGEE MOD 2**35 /JTD THEN CON 
C VEET THIS TO BLOATING POITID PSEUDORANDOM IlTO.rHBE UNIPOEEALY 

C DISTRIBUTED ON (0,1) . 

C PIRST STATEIffiNT PRODUCES EBSULT MOD 2**35 SIICE THE WATPOE 

C PORTRiiN IV SYSTEM ON THE IBM 7044 IGNOEBS OVERPLOW III INTC€£R 

C 1'IULTIPLICATION 

1 EN = ENl * AIPA 
ESTN = EN 

RSTN = RSTN *2.0«-*(-35) 
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0 THE EOLIOWliia PHOG-ruiJI 13 USED TO GOWERT iHaTEOM iriJl-IBERS 
UHIEOxHIIY DISTRIBUTED III THE IHTERYli (O.l) ^'0 EDTilSRS 
HRVIHG HOiaLii DISTRIBUTIOiT. 

0 IE (G-1 ),7,6,7 
7 H = H+1 
G = G+1 ^ 

RlTl = EH 
ESTHI = ESTH 
GO TO 1 

6 2 = SQET (-2.0*IiiOG(ESTHl))*SIGH;i 
Slfl = Z* 0 OS ( 6 ,2 83 18*ESTH ) + t!U 
3H = Z*SIF(6.28318*ESTH) + I;IU 
G = 0 
5 IT = H+1 
EHl = EF 
EETuni 
EHD. 
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